N82- 32923 

On cl as 
34555 

WA\S/V 

Technical Memorandum 83830 


AN EFFICIENT ROUTINE FOR 
INFRARED RADIATIVE TRANSFER 
IN A CLOUDY ATMOSPHERE 


""OSPHrPE (KASI, 5P p 

CSCL OUA 


G3/46 


Ming-Dah Chou and Louis Kouvarls 


OCTOBER 1981 



National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 


TM 83830 


AN EFFICIENT ROUTINE FOR INFRARED RADIATIVE 
TRANSFER IN A CLOUDY ATMOSPHERE 

by 

Ming-Dah Chou 

Goddard Laboratory for Atmospheric Sciences 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

and 

Louis Kouvaris 

Applied Research Corporation 
handover, Maryland 20785 
and 

Goddard Laboratory for Atmospheric Sciences 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 


October 1981 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 


“PAGE MISSING FROM AVAILABLE VERSION” 


PRECEDING PAGE BLANK NOT FILMED 


AN EFFICIENT ROUTINE FOR INFRARED RADIATIVE 

TRANSFER IN A CLOUDY ATMOSPHERE 

by 

Ming-Dah Chou 

Goddard Laboratory for Atmospheric Sciences 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

and 

Louis Kouvaris 

Applied Research Corporation 
handover, Maryland 20785 
and 

Goddard Laboratory for Atmospheric Sciences 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 

This report is the documentation of a FORTRAN program that calculates the atmospheric 
cooling rate and infrared fluxes for partly cloudy atmospheres based upon the fast but accurate 
methods of Chou and Arking (1978, 1980). The IR fluxes in the water bands and the 9.6 and 
15 /am bands are calculated at 15 levels ranging from 1.39mb to the surface. 

The program is generalized to accept any arbitrary atmospheric temperature and humidity 
profiles and clouds as input and return the cooling rate and fluxes as output. Sample calculations 
for various atmospheric profiles and cloud situations are demonstrated in the report. 

A digital magnetic tape containing the computer codes and the precomputed transmission 
functions of the radiation routine can be made available upon request. 


AN EFFICIENT ROUTINE FOR INFRARED RADIATIVE 


TRANSFER IN A CLOUDY ATMOSPHERE 


1 . INTRODUCTION 

Radiative heating and cooling play a crucial role on our climate system. The need for a fas 
yet accurate radiation routine has become more and more demanding as large-scale numerical 
models are now commonly used as a powerful tool to understand atmospheric processes and to 
study climatic problems. It has been a fact that the calculation of the radiative terms takes a 
large percentage, sometimes up to 90%, of the computing time spent in a numerical experiment. 
Usually, the percentage increases with the accuracy of a radiation routine that leads to the pain- 
ful choice of a crude radiation routine in long-term climate studies. 

Fast yet accurate radiation routines which include both water vapor and carbon dioxide bands 
have recently been developed by Chou and Arking (1978,^ 1980^). These routines for clear atmos- 
pheres are extended to include partly cloudy situations under the assumption that clouds are either 
black or gray. The computer codes of the radiation routine are described in detail in this report. 
Radiative transfer equations for fluxes are first presented in Section 2, simplifications of the 
transmission functions are then given in Section 3. Section 4 depicts the model structure and the 
approach to the vertical integration of the transfer equations. Section 5 describes the computer 
codings of the radiation routine. The computer program is listed in Appendix A with arrays and 
variables tabulated in Appendix B. Some sample calculations are demonstrated in Appendix C. 

2. RADIATIVE TRANSFER EQUATIONS 

In this radiation routine, the infrared spectrum is divided into four regions: the 15/im CO 2 
band, the 9.6jum O 3 band, the water vapor band center region and the band wing region. Clouds 

^Chou, M. D. and A. Arking, 1978. An infrared radiation routine for use in numerical atmospheric models. Pre- 
print, Third Conference on Atmospheric Radiation, American Meteorologies Society, 303-305. 

^Chou, M. D. and A. Arking, 1980: Computation of infrared cooling rates in the water vapor bands. J. Atmos- 
pheric Sciences, 37 , 855-867. 
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are treated as either black or gray, and scattering due to molecules and particulates is neglected. 
Fluxes in each spectral region are computed for clear and cloudy atmospheres. These fluxes are 
then weighted by the fraction of sky cover to obtain the total flux. 


For a clear atmosphere, the upward and downward fluxes at pressure level p integrated over 
the spectral range Av can be computed respectively from 


Fclr t (p) = 


= f dv 


B^(Ts) r^(p, Ps) + 


/' 

✓ D, 


dT^(p, p') 

B^(T(p')) — — dp 
dp' 


r 


(la) 


and 


C V dr^(p') I 

Fclr|(p)= j B/T(p')) dp' J (lb) 

where B^(T) is the blackbody flux at temperature T and wavenumber v and is equal to tt times 
the Planck radiance, T^(p, p') the diffuse transmittance between the levels p and p', and the sub- 
scripts s and t denote the surface and the top of the atmosphere, respectively. 


Because of the rapid variation of transmittance with wavenumber, the use of Eqs. (la) and 
(lb) to compute fluxes is a very time-consuming process. Simplifications of Eqs. (la) and (l,b) 
are applied to different spectral bands. In the 9.6 and ISjtim bands, the transmittance Tp is re- 
placed by its average value over which reduces the equations for fluxes to 


Fch t (p) = B(Ts) r(p, p^) + f 

•'Ps 


P dr(p, p') 

B(T(p')) p- dp' 

dp 


and 


/*P dr(p, p') 

Fclr |(p) = / B(T(p')) dp', 

•'Pt 


dp' 


where 




B(T) = f B^(T)di;. 
Jt^ 


(2a) 


(2b) 


(2c) 


2 


and 


r(p, p') = 






T^(p, p') di^. 


( 2 d) 


In the water vapor bands, the flux computations are simplified by integrating Eqs. (la) and 
(lb) by parts and reducing them to 


pclr t (p) - B(T(p)) + G(p, p,, T.) - G(p, Ps, T(pJ) + 




Ps> HPsJ) + I 3G(p, p', T(p'))/0T dT(p'), (3a) 

^(P) 


and 


Fc1t4(p) = B(T(p)) - G(p, pt, T(p 


/•T(P,) 

^t)) + / ^ 

•'t(p) 


3G(p, p', T(p'))/3T dT(p'), 


(3b) 


where 


and 


G(p, pM) = f r^(p, p') B^(T) di.. 
J/Sp 


3G(p, p', T)/3T = f r^(p, p') 3Bj,(T)/3T di^. 
J Av 


(3 c) 


(3d) 


It is noticed that the Planck function varies significantly within a wide spectral interval. The 
use of Eqs. (2a) and (2b) to compute fluxes would introduce significant errors unless the water 
vapor bands are divided into a number of intervals (^ 10). With the G-function computed from 
Eq. (3c) using the line-by-line method, the effect of the variation of Planck function with wave- 
number is correctly taken into account. The reason that Eqs. (3a) and (3b) are applied only to 
the water vapor bands is that we have not yet developed methods for precomputing the G- 
function in the CO 2 and O 3 bands. This will become clear in Section 3. 


For a cloudy atmosphere, Eq. (2a) is used to compute fluxes in the 9.6 and 15)um bands except 
that the subscript s is replaced by CT for upward fluxes and CB for downward fluxes, where CT de- 
notes the cloud top and CB the cloud base. Similarly, Eq. (3a) is used to compute fluxes in the water 
vapor bands with the subscript s replaced by CT and CB respectively for upward and downward fluxes. 
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For a partly cloudy atmosphere, the infrared flux is computed from 

n 

F(p) = (1 -C)Fdr(p)+^CiFiCld(p), (4) 

1 = 1 

where F^^d jg the flux for a cloudy atmosphere, C the total fractional cloud cover, Cj the cloud 
cover of type i, and n the number of cloud types. The cooling rate is proportional to the flux 
divergence given by 

dT g d 

- — = (f; - Ft), (5) 

dt Cp dp 

where Cp is the heat capacity of air and g the acceleration of gravity. 

3. COMPUTATION OF TRANSMISSION FUNCTIONS 

The most difficult part in computing radiative fluxes in the atmosphere is the integration 
over the entire spectrum. Usually the absorption coefficient changes several orders of magnitude 
within a spectral interval of 100 cm"*, and the transmittance in Eq. (1) must be computed at 
every ^0.01 cm"* interval, which results in a total of 10^-10® points along the spectrum. Since 
we are interested in the total flux but not the flux at individual wavenumber, the point-by-point 
calculation is commonly replaced by wide band (or emissivity) or narrow band approximations. 

In this section, we describe the method used in our radiation routine to eliminate the difficulty 
in computing transmission functions in various parts of the spectrum. 

a. Basic Equations 

The beam transmittance between any two pressure levels pj and P 2 at wavenumber v in the 
direction from the vertical 6 = cos"*/a is 

= exp ^ Uj,(i)/M^ (6) 

where u^(i) is the optical depth of the ith absorber given by 


4 


k^(i) q(i)dp, 


(7) 


1 /"P2 

u^(i) = - / 

g 


kj^ the absorption coefficient and q the mixing ratio. 


The diffuse transmittance in Eqs. (1) and (3c) is related to the beam transmittance by 


= 2 


J Ty(ju) p dn, 
0 


( 8 ) 


b. Absorption Coefficients 


The absorbers considered in this radiation routine include water vapor, carbon dioxide, and 
ozone. The parameters needed for computing the molecular line absorption are the position, Vq, 
strength, S, width, a, and shape of molecular lines. Using the parameters compiled by McClatchey 
et al. (1973^) and assuming a Voigt line function, the absorption coefficients are computed using 
the line-by-line method at every 0.01 cm*' interval for individual absorbers. 


The e-type continuum absorption due to water vapor in the spectral region between 340 and 
1200 cm*' is also included. The regression curve given in Roberts et al. (1976^) which fits labor- 
atory data is used to compute the absorption coefficient at 296 K. 

c. Scaling Approximation for Water Vapor Absorption 

In the water vapor bands, the efficiency of the flux computation can be further increased 
by taking advantage of the nature of water vapor distribution. Generally the water vapor de- 
creases logarithmically with height. The lower troposphere is too opaque to have significant cool- 
ing for the spectral intervals near the center of the absorption bands, whereas the upper tropos- 
phere is too transparent to have significant cooling for the band wings. As shown in Chou and 
Arking (1980) that the cooling rate has a relatively narrow vertical profile that peaks at the upper 

^McClatchey, R. A., W. S. Benedict, S. A. Clough, D. E. Burch, R. F. Calfee, K. Fox, L. S. Rothman and J. S. Caring, 
1973: AFCRL atmospheric absorption line parameters compilation. Environ. Res. Pap., No. 434, AFCRL-TR-73- 
0096. 

^Roberts, R. E., J. E. A. Selby, and L. M. Biberman, 1976: Infrared continuum absorption by atmospheric water 
vapor in the 8-12jum window. Applied Optics, 15, 2085-2090. 
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troposphere for the band center region and at the lower troposphere for the band wing region. 

By scaling the absorption coefficient (or alternatively the water vapor amount) separately for the 
band center region and the band wing region based on the temperature and pressure at the height 
where the cooling rate is a maximum, the difficulty in computing the absorption coefficient for 
various temperature and pressure can be greatly reduced. 


At the levels where cooling is significant, absorption quickly saturates near the center of most 
of the molecular lines, and the cooling rate is mainly due to the region m line wings. In the far 
wings, where I n - I > a, the water vapor molecular line absorption, 2^, can be written as 


£^(p, T) = £^(p,, T,) — R(T, T,) 

Pr 


where R is the spectrally averaged value of R^ given by 


R^(T, T, 


= ( 


r Sj(T) Oj (p„ T,)- | 

i L -I 


T/ 


pidr) (Pp Tj)" 

4-L 0^ - - 


and the index i denotes absorption lines. 


(9) 


( 10 ) 


In Eq. (9), the absorption coefficient at the reference conditions, 2j,(pj, T^), is “exactly” 
computed using line-by-line method. The error in 2^(p, T) using the far-wing scaling approxima- 
tion depends upon the departure of p and T from pj. and T^. Since the water vapor cooling pro- 
file spans only a narrow region in the vertical, the error in cooling rate arising from the use of 
Eq. (9) is minimal if we choose the reference temperature and pressure to be in the region of 
maximum cooling. Table 1 shows the spectral regions, the reference temperature and pressure, 
and the mean values of R^ for the band center region and the band wing region. The reference 
temperature and pressure are chosen to correspond to conditions at the level of maximum cooling 
rate in the U.S. Standard Atmosphere. In Table 1 , the temperature scahng term R is given only 
at two temperatures. Values of R at other temperatures are computed using a quadratic fit to the 
R-values at T = Tj. ± 40K and T = Tj. (for which R = 1). 
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Table 1 


The reference temperature T^ and pressure Pj, are specified for the center and the 
wing regions of the water vapor absorption bands. The effect of temperature 
on the absorption coefficient is indicated by R(T, T^), which is the mean 
value of the function R^(T,Tf) defined by Eq. (10). 



Wavenumber (cm"') 

Tr(K) 

Pr(mb) 

R(Tr -40. Tr) 

R(T, +40, Tr) 

Center 

0-340 

1380-1900 

225 

275 

0.90 

1.16 

Wing 

340-580 

760-980 

1100-1380 

1900-3000 

256 

550 

0.58 

1.78 


d. Diffuse Transmittance in the 15/nm Band 

The absorption in the 15 /am band from 580 to 760cm"' is due primarily to CO 2 and second- 
arily to water vapor. The overlapping of diffuse transmittances is approximated by 

’"is/a “ '’'C02 ■ ’’e ■ ’’e’ ( ^ 

where Tqq^, rj, and are the transmittances averaged over the 15/am band due to CO 2 molecular 

line absorption, H 2 O line absorption and H 2 O e-type absorption, respectively. 

The atmospheric CO 2 content can practically be considered as constant, and the transmit- 
tance TCO 2 between any pressure levels is a function of temperature only. Since the relative 
change in atmospheric temperature is not large, Arking (1976^) used a linear expansion method to 
correct for small variation in the beam transmittance for temperature retrievals. This method in- 
troduces an error ^0.005 in the mean diffuse transmittance averaged over the entire 15 pm band 
(Chou and Arking, 1978). In this infrared radiation routine, the transmittance Tqq^ between pj 

and p 2 is computed using the linear expansion method from 

m 

’■c02(Pi> P 2 ) = 'rc02(Pi’P2) H aj(pi, P 2 ) ATj, (12) 

i=l 

^Arking, A., 1976: The proper weighting function for deriving temperatures from satellite measured radiances. 

GSFC X-document X-91 1-76-48. 
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where the superscript o denotes the transmittance for a reference temperature profile, m the 
layers between pj and P 2 , and AT the deviation of temperature from the reference temperature. 
The coefficient a^ is given by 

ai(Pi.P2) = 9 t(P], P2)/3Tj. (13) 

The e-type continuum absorption coefficient, C^, increases with vapor pressure but decreases 
with temperature according to (see Roberts et al., 1976) 

r 1800 ISOOI 

C^(T) = C^(To) e • exp 1^- — + — J , (14) 

where C^(Tq) is the absorption coefficient at p = 1 atm and T = Tq, and e the vapor pressure in 
atms. 


With the use of Eq. (14) and the far-wing scaling approximation, Eq. (9), the transmittances 
Tg and Tg are computed from 


Tg(w) = 


Tg(u) 


lur 

Jo 

= 1 fa. r 

Jav Jo 


exp(-e^(Pr, Tr)w/ju) p dp. 


exp(-C^(T(j)u/p) p dp. 


•' 0 

where the scaled water vapor amounts are given by 

1 /*P2 / P 


w(Pi, P 2 ) 


1 rP2 / 

u(Pi, P 2 ) = — / e • exp ( - 

g Jpi ^ 


= - f (—) R(T, T,) q dp, 
g Jpi \Pr/ 

1800 1800 


T T 


q dp. 


(15) 


(16) 


(17) 


(18) 


Values of 1 ^ 2 )’ ^(Pi> P 2 )’ T’g(w), and Tg(u) are precomputed using the line-by-line 

method and stored in tables for later use. The reference conditions are chosen to be the same as 
the water vapor band wing region, i.e., Pj. = 550mb and Tj. = 256K. Values of R are approximated 
by that given in Table 1 for the band wing region. 
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e. Diffuse Transmittance in the 9.6um Band 


In the 9.6Mm band from 980 to 1 lOOcm"^ the absorption is mostly due to ozone in the 
stratosphere. The water vapor e-type absorption contributes somewhat to the cooling in the 
lower tropical atmosphere, but the effect of line absorption is negligible. Applying the multipli- 
cation approximation, the diffuse transmittance is computed from 

T9.6M ~ ^03 ' ^e- 

The ozone transmittance, Tq^, is precomputed using the line-by-line method for the three 
latitude zones, 0-30, 30-60, 60-90°N. Annual mean ozone content and temperature profile at 
15, 45, and 75°N are used to compute for the respective latitude zones. No attempts have 

been made to correct for the change in due to changes in the ozone content and temperature 

profile. The transmittance Tq is precomputed from Eq. (16) and stored in a table as a function 
of the scaled water vapor amount, u. 

f. The G-Function for the Water Vapor Bands 

Utilizing the far-wing scaling approximation, the water vapor transmittance between pj and 
P 2 can be written as 

T,;(Pi, P 2 ) - ^ f exp - [(S^(Pp Tf) w(pi, P 2 ) + Cy(To) u(pi, P2 ))/m1m dja; (20) 

Jo 

where w(pj, P 2 ) and u(pj, P 2 ) are given by Eqs. (17) and (18). The error introduced in T^(Pi, P 2 ) 
by using Eq. (20) was shown in Chou and Arking (1980) to be small with a magnitude :^0.01. 

Since 2^(Pp T,.) and C^(Tp) are functions of wavenumbers only, the G-function defined in 
Eq. (3c) depends on w, u, and T, but not on the detailed temperature and humidity profiles. 

The function G(p, p', T) in Eqs. (3a) and (3b) can then be approximated by G(w, u, T) with w 
and u being the scaled water vapor amounts in the column between p and p'. It is now possible 
to precompute the three-dimensional function of G(w, u, T) using the line-by-line method and 
the results stored in tables. The G-function defined in Eq. (3c) is computed for the band wing 
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region and the band center region with spectral interval Ai' listed in Table 1 . A table look-up for 
the G-function makes the computations of fluxes very fast from Eqs. (3a) and (3b). 

4. VERTICAL INTEGRATION FOR FLUXES 

a. Model Structure 

Since the functions B, t, and G in Eqs. (2) and (3) are highly non-linear with respect to 
height, the vertical integrations cannot be performed analytically. The continuous atmosphere is 
therefore discretized by dividing it into a number of layers, and the integrations approximated by 
summations. The fluxes are computed at 15 levels (designated as flux levels hereafter) separated 
by A In p = 0.828 above the 200 mb level and Ap = 100 mb below. Since the layers are too thick 
(in terms of the differences in the Planck and transmission functions at the boundaries) for flux 
computations, each layer is further divided into 3 sublayers with equal A In p above the 200 mb 
level and equal Ap below. Thus, there are 43 sublevels and 42 sublayers in the vertical. Fig. 1 
illustrates the pressure at these levels. 

b. Temperature and Humidity Interpolations 

Given temperature and humidity at any pressure levels as inputs to the radiation routine, 
they are interpolated (or extrapolated) to the 43 sublevels and the middle of each of the 42 sub- 
layers. Temperature and the logarithm of specific humidity are interpolated linearly in In p. 

c. Vertical Integrations 

1 

Let j be the index for sublevels increasing downward and j + -j for the sublayer immediately 
below, then the following approximations are made in computing fluxes at the level pj in the 9.6 
and 15jum bands (see Fig. 1 for the locations of i and j), 

B(T(p')) dr(pj, p') ^ B(Tj^) [r(pj, Pj.i) -T(pj, pj)] , j = i + 1, . . . , 43, for upward flux; (21a) 
~ B(Tj+^) [t(Pp Pj+i) -t(Pj, pj)], j = 1, . . . , i - 1, for downward flux. (21b) 

Similarly, the following approximations are made in computing water vapor infrared fluxes. 
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PRESSURE FLUX LEVEL SUBLEVEL 

(mb) IP IX 
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- 24 



Figure 1. Model Vertical Structure. Sublevels are Used for 
Computing Fluxes at the Flux Levels. 


11 


ORIGSr^AL PAGE SS 

OF POOR QUALITY 

[3G(Pi,p',T(p'))/3T(p')] dT(p') 

~ G(w(pj, pj^), u(pi, pj^), Tj) 

- G(w(pj, pjj^), u(pj, pj^), Tj.j), j = i + 1, . . . , 43, 

for upward flux; (22a) 

~ G(w(pj+,^, Pi), u(pj+,^, Pi), Tj) 

- G(w(pj+,^, pj), u(pj+v 4 . Pi), Tj+i),j = 1, . . . , i - 1, 

for downward flux. (22b) 

5. INFRARED RADIATION CODE 

The program is divided into 6 sections. The main routine is for specifying the input data and 
for returning the fluxes and the cooling rate . The fixed data that accompany the program are read 
in the subroutine INDATA. The third routine GAMOTO interpolates the input data of cloud heights, 
temperature, and humidity to the model levels and computes the scaled water vapor amounts for 
each layer. The subroutine IRH20 interpolates the G-function in the water vapor bands which is used 
in the subroutine FDIVER to compute fluxes in the water vapor bands. Finally, the subroutine 
RC0203 computes the fluxes in the 9.6 and 1 5/zm bands. Subroutines IRH20 and RC0203 are in- 
dependent of each other. In the following, variables and arrays in the routines are all denoted by 
capital letters. 

a. Main Program 

The main program is to accept the input data and return the computed fluxes and cooling rate. 
The users must specify the following input data 


Input Data 

Unit 

Co 

Pressure 

mb 

PA 

Temperature 

K 

TA 

Humidity 

gm/gm 

WA 

Surface temperature 

K 

TS 

Surface pressure 

mb 

PS 
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Input Data 

Unit 

Code 

Cloud top pressure 

mb 

CLTOP 

Cloud base pressure 

mb 

CLBOT 

Fractional cloud amount 


CSS 

Number of pressure levels 

NDATA 

Number of cloud types 


NC 

Latitude 

degree 

RLAT 

The user’s input of pressure, temperature and humidity define the atmospheric profiles. These 

data must be arranged in order of increasing pressure. The surface pressure, PS, must be greater than 

900mb. Different cloud types are 

treated as non-overlapping. 

For cloud-free atmospheres, the 

parameter NC is set to zero. 



The following parameters are 

the outputs of the infrared 

radiation routine: 

Parameter 

Unit 

Code 

Downward flux in the 
water vapor bands 

ergs/(cm^ sec) 

FLUXHD 

Upward flux in the 
water vapor bands 

ergs/(cm^ sec) 

FLUXHU 

Downward flux in the 
9.6 and 15/Ltm bands 

ergs/(cm^ sec) 

FLUXCD 

Upward flux in the 
9.6 and 1 5/xm bands 

eigs/(cm2 sec) 

FLUXCU 

Cooling rate 

°C/day 

COOLR 


Fluxes are computed at the 15 levels denoted by the running variable IP in Fig. 1, and cool- 
ing rates are computed for the layers defined by the flux levels. The spectral ranges for water 
vapor bands are listed in Table 1. The cooling rate is computed from Eq. (5) with Cp = 

1.003 X 10’ ergs/(gm K). 
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b. Subroutine INDATA 


This subroutine reads in the model pressure sublevels and the precomputed data for trans- 
mission and G-functions. As mentioned in Sec. 4, fluxes are computed at 15 levels (NP = 15) 
and cooling rates for 14 layers. For proper vertical integration of the transfer equations, each 
layer is further divided into 3 sublayers (NSB = 3). Therefore we have 43 sublevels (NPT = 43) 
defining 42 sublayers. The pressure of the first 40 sublevels (PUI) are read in this subroutine 
and fixed in the model. The last 3 sublevels are computed in subroutine GAMOTO, which de- 
pends on the surface pressure. 

The G-function is stored in the 3-dimensional array GFH20. The first index runs from 1 to 
25 (NT = 25) corresponding to temperatures of 190K (TEMPI = 190) to 31 OK incremented in 
steps of 5K (DT = 5). The second index corresponding to log w runs from 1 to 22 (NW = 22) 
where the first value corresponds to log w = -6 (SW = -6) for the band center absorption and to 
~5.4 (SW = - 5.4) for the band wing absorption. It is incremented in steps of A log w = 0.3 (DW = 
0.3). The third index corresponding to u runs from 1 to 11 (NUPl = 11). Index values from 1 
to 10 refer to the band wing region given at Au = 0.006 (DU = 0.006) intervals starting from 
u = 0. The last index value is for the band center region which does not include e-type absorption. 

The following table summarizes the information inGFH20 (IT, IW, lU) for computing fluxes 
in the water vapor bands. The units are K for T and gm/cm^ for w and u. 

Table 2 


Information in the Precomputed G-Function 


Index 

Corresponding Parameter 

Initial Value 

Increment 

Band-Center 

Band-Wing 

IT 

T 

190 

190 

5 

IW 

i 

log w 

-6.0 

-5.4 

0.3 

lU 

u 

- 

0 

0.006 


14 














In the 15ptm band, the CO 2 transmission function, Tqq^, is stored in TRE as a two- 
dimensional array. The first index refers to sublevels and the second index to flux levels. The 
coefficients for the linear expansion, aj in Eq. (12), are stored in the one-dimensional array FSN 
only for the region above the 200 mb level. FSN is ordered in a way as shown in Fig. 2. Values 
of TRE and FSN are computed for the U.S. Standard Atmosphere. The transmittance due to 
water vapor line absorption, rg, is stored in TAUW for 65 values of log w (NW = 65) starting 
from -4.0 (SW = -4) incremented by 0.1 (DW = 0.1). The transmittance due to e-type ab- 
sorption, Tg, is stored in TAUU. Values of TAUU are given at A log u = 0.04 (DU = 0.04) inter- 
vals starting from log u = -3.301 (SU = -3.301). 

For the 9.6)um band, the ozone transmittance, Tq^ is given by the three-dimensional array, 
TOZON. The first index refers to sublevels, the second index to flux levels, and the last to 


PRESSURE 

(mb) 

1.39 


1.83 
2 41 


i FSN I FSN 
1(192) 1(191) 


FSN(193) 


115.1 

151.8 

200 


i FSN 
1(1140) 


i FSN(1137) 


FSN(1139) 


SUBLEVEL 


FSN(2) 


FSN(I) 


FSN(3) 


4 FSN(189) 


FSN(190) 


3 

4 


19 


Figure 2. Ordering of the Coefficients FSN’s for the Correction of CO 2 
Transmittance in the Stratosphere. 
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latitude zones. In computing ozone transmittance, annual mean ozone mixing ratios representative 
of 30°-latitude zones are used. The last index (LAT) is equal to 1,2, and 3 for the latitude 
zones 0-30, 30-60, and 60-90°, respectively. The transmittance due to e-type absorption, Tg, is 
stored in TWIN as a function of log u. It is given at A log u = 0.04 (DU = 0.04) intervals start- 
ing from log u = -3.301 (SU = -3.301). 

The following table summarizes the information in TAUW, TAUU, and TWIN. The units for 
w and u are gm/cm^. 


Table 3 


Information in the Precomputed Transmittances due to Water Vapor Absorption 


Transmittance 

Spectral Band 

Index 

Corresponding Parameter 

Initial Value 

Increment 

TAUW 

15/im 

IW 

log w 

-4.0 

0.1 

TAUU 

1 5jum 

lU 

log u 

-3.301 

0.04 

TWIN 

9.6jum 

lU 

log u 

-3.301 

0.04 


c. Subroutine GAMOTO 

This subroutine is to 

i) Interpolate user’s atmospheric profiles onto the model grids; 

ii) Compute scaled water vapor amounts; 

iii) Compute indices for cloud top and cloud base; 

iv) Call subroutines IRH20 and RC0203. 

The two sublevels between 900 mb and the surface are computed by dividing the layer into 
three sublayers with equal Ap, i.e. Ap = (p^ - 900)/3. The center pressure (PAI) and the air mass 
(DPI) of each sublayer are then computed. The temperature at each sublevel (TUI) and the tem- 
perature and humidity in the middle of each sublayer (TAI and RAWI) are interpolated from. the 
input values of pressure (PA), temperature (TA), and specific humidity (WA). Temperature and 
the logarithm of the specific humidity are interpolated linearly in In p. The same parameter 
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RAWI designated for specific humidity in line 01800 is later redesignated for water vapor content 
in line 01840. 

The scaled water vapor amounts w (WEAR) and u (UBAR) in a column above any pressure 
level are computed in lines 01910-02000 from the above parameters using Eqs. (17) and (18). 

The reference pressure p^ is coded as PRE and is equal to 225 mb for the band center region and 
to 550 mb for the band wing region. The mean effect of temperature on the absorption coeffi- 
cient, R, coded as SZ2, is computed from the R-values in Table 1 using a quadratic fit. The 
coefficients for the quadratic fit are coded as RR. The index KB for the arrays WEAR, RR, and 
PRE is equal to 1 for the band center absorption and equal to 2 for the band wing absorption. 

In computing values of u, the following relationship is used 

eq/g dp = ( — (1.02 q) dp' = (1.02/630) p'q2 dp'. 

\0.622 P/ 

where the units are dynes/cm^ for p and mb for p', and P is the standard pressure equal to 
1013mb. The temperature Tq in Eq. (18) is equal to 296 K. 

The pressure corresponding to the cloud levels are changed to the nearest sublevel numbers... 
of the model grids. The cloud top pressure (CLTOP) is changed to NXT, the cloud base (CLBOT) 
to NXB, and the cloud amount (CSS) is just redefined to CS. 

d. Subroutine 1RH20 

This subroutine computes the terms on the right hand sides of Eqs. (3a) and (3b), which are 
then used in the subroutine FDIVER for computing fluxes in the water vapor bands. The Planck 
flux is either coded as SH at each sublevels or as FLUXUS at the surface. G(p, p^, T(pj)) is 
coded as SHT, [G(p, Pg, T^) - G(p, p^, Tfpj))] as SHS, and (8G/3T) dT as SG. These terms are 
all computed from the discretized G-values using linear interpolation. 

Before discussing the computation of the Planck flux, it is noticed from Eq. (3c) that 
B(T) = G(w = 0, u = 0, T) = G(wi,Ui,T) + 6(T) 
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where Wj and Uj are, respectively, the starting values of w and u for the precomputed G-function, 
and 5 (DELTA) the correction term. 

The subroutine starts in lines 02560-02670 by interpolating the G-values, corresponding to 
the surface temperature and the values of Wj and Uj , for the band wing region (NZ = 1) and 
the band center region (NZ = NUPl). The upward flux at the surface (FLUXUS) is then the sum 
of the two G-values and the term DELTA. A similar approach is later applied in computing the 
Planck flux at each sublevel (SH) in lines 03190-03260. 

Following the computation of FLUXUS is the computations of the index value for the tem- 
perature nearest to the sublevel temperature and the difference between the sublevel temperature 
and the temperature at the nearest index value. They are stored in IH and DH for later use. 

The rest of the subroutine is used to compute SG with SHT and SHS as by-products. Given 
values of w, u, and T, the G-function is computed from 

3G 3G 3G 

G(w, u, T) = G(Wj, Uj, Tj^) + (log w - log wj) + — (u - m) + — (T - T^) (23) 

^ 9 log w 9u J 9T 

where i (KW), j (KU), and k (IT) are the indices for the precomputed G-function nearest to log w, u, 
and T, respectively. The term log w - log Wj is coded as FW, u - Uj as FU, and T - Tj^ as DH. The 
computation of G-function using (23) is coded in lines 03030-03080. Values of SG are then com- 
puted from Eq. (22) and that of SHS from 

9G 

G(w, u, Tj) - G(w, u, T(ps)) = — (Tj - Tfp^)) 

3T 

~ [G(Wj, Uj, Tj^+i) -G(wj, uj, Tj(.)]/(Tk+i - T^) • (T^ -T(pg)). 

Values of SG are computed in lines 03080 and 03120, SHT in 03090 and SHS in 03100. 

e. Subroutine FDIVER 

This subroutine computes upward and downward fluxes in the water vapor bands (FLUXU 
and FLUXD) from Eqs. (3a) and (3b) using the parameters SH, FLUXUS, SHT, SHS, and SG 
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computed in the subroutine IRH20. Treatments of various cloud situations are given in Sec. 
4g- 

f. Subroutine RC0203 

This subroutine computes the upward and downward fluxes in the 9.6pm band (FLUXUW 
and FLUXDW), the 15pm band (FLUXUC and FLUXDC), and the total in these two absorption 
bands (FLUXU and FLUXD). 

The Planck fluxes integrated over the 9.6pm band and the 15pm band are interpolated in 
lines 04340-04470 from the precomputed Planck fluxes. At the middle of each sublayer, the 
Planck fluxes are coded as BAI for the 15pm band and BWl for the 9.6pm band. The corre- 
sponding Planck fluxes at the surface are coded as BS and BSS, respectively. The precomputed 
Planck flux in the 9.6pm band is stored in BLKWIN and that in the 15pm band in BLKC02 at 
5K intervals starting from 190K. 

The transmittances due to the ozone line absorption (TOZON) and the CO 2 absorption 
(TRE) are first transferred to the arrays TXUF and TAUF, respectively, in lines 04520-04550. 
The correction for the CO 2 transmittance above the 200 mb level due to the departure of the 
temperature profile from the U.S. Standard Atmosphere (TR) is then coded in lines 04570- 
04780. (Note that no correction is made if the parameter LINEXP is set to FALSE). The 
transmittance due to water vapor line absorption (X) is computed in lines 04830-04960, and 
the transmittances due to e-type absorption are computed in lines 04980-05090 for both the 
15pm band (Y) and the 9.6pm band (YY). Finally, the total transmittances are computed from 
Eqs. (11) and (19) and are coded in lines 05110 and 05120. 

With the computed Planck and transmission functions, fluxes in the 9.6 and 15pm bands 
are computed from Eqs. (2a) and (2b) in the rest of the subroutine. The treatments of various 
clouds are explained below. 
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g. Treatments of Clouds 


Clouds are assumed to be black in the infrared. The routine can be easily extended to a 
gray cloud with a constant emissivity. It can be shown that a gray cloud with emissivity e and 
amount C is radiatively equivalent to a black cloud with amount eC. Although the radiative 
transfer equations used for computing fluxes in the water vapor bands and in the 9.6 and 15/im 
bands are different, the treatments of clouds are similar. Fig. 3 illustrates an atmosphere with 5 
possible cloud situations. The indices NXT and NXB refer to the cloud top and cloud base, re- 
spectively, and C the fractional cloud cover. In the H 2 O bands, the computation of the down- 
ward flux at the pressure level IP is coded in lines 03450-03710 according to the following 
relationships, 

IS-l 

FLUXD(IP) = SH(IP) - SHT(IP) + SG(I, IP) 

1 = 1 

for clear atmospheres or IS < NXT 
IS-l 

= SH(IP) + 51 SG(I, IP) for IS > NXB 

I=NXB 

= SH(IP) for IS = NXB or NXB < IS < NXT 

The total flux is the sum of these fluxes weighted by respective fractions of sky cover. For the 
case shown in Fig. 3, the fraction of the clear sky for the downward flux at IP is 1 - (C 2 + C 3 + 
C 4 ) which is coded as 1 - SUMC. 

At each flux level, IP, a one-dimensional array of FLUX is computed from 

IS-l 

FLUX(IX) = SG(I, IP), IX = 1 IS - 1 , 

I = IX 

which is used for computing fluxes in various cloud situations. Computations of upward fluxes 
are treated similarly and are not shown here. 
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SUBLEVEL 



id Base, Respectively, and the C’s are Cloud Amount. 







With the aid of Fig. 3, the computation of the upward flux in the 15pm band coded in 


lines 05160-05560 is explained below. The upward flux at the level IP is computed from 

FLUXUC(IP) = BS * TAUF(NPT, IP) 

NPT 

+ * [TAUFd - 1, IP) - TAUFd, IP)] 

I=IS+1 

for a clear atmosphere or IS > NXB 
= BAKNXT + 1) * TAUF(NXT, IP) 

NXT 

+ ^ BAI(I) * [TAUFd - 1, IP) - TAUFd, IP) 

I=IS+1 

for IS < NXT 

= BAKIS +1) for IS = NXT or NXB < IS < NXT 
An intermediate one-dimensional array FLUXC is computed from 

IX 

FLUXC(IX) = BAI(I) * [TAUFd - 1, IP) - TAUFd, IP)] , 

I=IS+1 

IX = IS + 1, ... , NPT 

and is used for computing fluxes in different cloud situations. The fraction of clear sky for the 
case shown in Fig. 3 is 1 - (Cj + C4 + C5) and is also coded as 1 - SUMC. 
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APPENDIX A 


PROGRAM LIST 


DIMENSION FLUXHUa5>,FLUXHDU5;,FUJXCUa5>,FLUXCDa5/ 

DIMENSION FXWAVU5j,FX15Ca5;,COOLRa5;,PUU43j 
COMMON/FLUXH/ FLUXHU , FLUXHD 
COMMON/FLUXC/FLUXCJ , FLUXCD 

COMMON/AMHN/TA(.55> ,PAc55j ,WA^55) .CSSUOj , NDATA , RLAT , TS , PS 
COMMON/PRESS/PUI 

COMMON /INIT/NUP1,NP,NSB,NC,NPT,NW,NT 
COMMON / CLOU / CLTOP U 0 CLBOT U 0 ; 

READ (. 4 , 2 6 ; NDATA , NC , RLAT , TS , PS 
26 FORMATa2,2X,I2,2X,F4.2,2X,F5.1,2X,F6.1> 

READ I 4, 32 J DUMMY 

READ(,4 , 21 nPAU > , 1=1 .NDATA; 

21 FORMAT15F6.1; 

READ 14, 32; DUMMY 
READi4,23mTAU; ,1=1, NDATA; 

23 FORMAT18F7.2; 

READ 1 4, 32} DUMMY 
READl4,24;iWAU; ,1=1, NDATA; 

24 FORMAT15E9.7; 

READ 14,32; DUMMY 

READi4,25;iCLTOPiI; .CLBOTiI; ,CSS U; , 1=1 , 5; 

25 FORMATl5l2F4.0,F4.2;; 

PRINT 10 

10 FORMAT 112X,' NDATA', 2X,'NC' ,2X,'RLAT' ,3X, ' TS ' , 6X, 'PS ' ; 

PRINT 12,NDATA,NC,RLAT,TS,PS 

12 FORMATll3X,I2,4X,I2,2X,F4. 1,2X,F5. 1,2X,F6.1; 

PRINT 14 

14 FORMAT l/,l OX, 'PRES SURE LEVELSlMB; OF SAMPLE PROFILE'; 

PRINT 15, iPAlI; ,1=1, NDATA; 

15 FORMATl5l2X,F6.1,6X;; 

PRINT 16 

16 F0RMATl/,10X, 'TEMPERATURES 1K;AT THE ABOVE PRESSURE LEVELS'; 

PRINT 11, ITAU; ,1=1, NDATA; 

11 FORMATl5llX,F8.2,5X;; 

PRINT 17 

17 FORMAT l/,10X,' HUMIDITY IGM/GM; AT THE ABOVE PRESSURE LEVELS'; 

PRINT 13, iWAlI; ,1=1, NDATA; 

13 FORMATl5llX,E9.3,4X;; 

PRINT 19 

19 FORMAT l/,l OX, 'CLOUD TOP ', 2X ,' CLOUD BOTTOM' ,2X,' CLOUD PERCENT'; 
PRINT 18,lCLTOPlIC;,CLbOTlIC; ,CSS lIC; , IC=1 , NC; 

18 FORMATll4X,F4.0,7X,F4.0,9X,F4.2; 

32 FORMAT 1A1,/,A1; 

CALL INDATA 
CALL GAMOTO 
PRINT 3 

3 FORMATl//, 3X, 'PRESSURE' ,4X, 'FLUXES IN H20 BANDS' ,3X, ' FLUXES IN 9. 
$6 & 15 MICRON BANDS'; 

PRINT 4 

4 FORMAT15X,' IMB; ' ,7X,' lERGS/CM**2/SEC;' , IIX, ' lERGS/CM**2/SEC; ' ; 
PRINT 5 

5 FORMAT ll5X,' DOWN' ,8X,'UP' , 14X, 'DOWN' , 8X,'UP' ; 

DO 2 IP=1,NP 

IPM1=IP-1 

IDD=lIPMl*NSB;+l 

2 PRINT l,PUIlIDD; ,FLUXHDlIP; ,FLUXHUlIP; ,FLUXCDlIP; ,FLUXCUlIP; 

1 FORMATllX,F8.2,2F12.2,4X,2F12.2; 

F,X=3. 6*24./ 10030. *.98 
DO 69 IP=1,NP 

FXW AV a P ; “ FLU XHD U P ; - FLU XHU U P ; 


IR 

OUOlO 

IR 

00020 

IR 

00030 

IR 

00040 

IR 

00050 

IR 

00060 

IR 

00070 

IR 

00080 

IR 

00090 

IR 

00100 

IR 

00110 

IR 

00120 

IR 

00130 

IR 

00140 

IR 

00150 

IR 

00160 

IR 

00170 

IR 

001 80 

IR 

00190 

IR 

00200 

IR 

00210 

IR 

00220 

IR 

00230 

IR 

00240 

IR 

00250 

IR 

00260 

IR 

00270 

IR 

00280 

IR 

00290 

IR 

00300 

IR 

00310 

IR 

00320 

IR 

00330 

IR 

00340 

IR 

00350 

IR 

00360 

IR 

00370 

IR 

00380 

IR 

00390 

IR 

00400 

IR 

00410 

IR 

00420 

IR 

00430 

IR 

00440 

IR 

00450 

IR 

00460 

IR 

00470 

IR 

00480 

IR 

00490 

IR 

00500 

IR 

00510 

IR 

00520 

IR 

00530 

IR 

00540 

IR 

00550 

IR 

00560 

IR 

0‘,/570 

IR 

00580 

IR 

00590 

IR 

00600 
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69 

FX15CaP>=FLUXCDUP>-FLUXCUaP> 

IR 

00610 


PRINT 6 

IR 

00620 

6 

FORiiATv//16X/PRESSURE' ,10;;/COOLING RATE' ; 

IR 

00630 


PRINT 7 

IR 

00640 

7 

F0RMATU6X / ^MBAR/ ,7X,' ^DEGREE CELCIUS/DAY/ ; 

IR 

00650 


PRINT 8 

IR 

00660 

8 

FORMATU 2X , ' FROM' , 6X/ TO' y 

IR 

0u670 


DO 90 IP=2,NP 

IR 

00680 


IPM1=IP-1 

IR 

00690 


IDD=UPM1*NSB;+1 

IR 

00700 


IDE=IDD-NSB 

IR 

00710 


DP«=PU I UDD > -PU 1 UDE y 

IR 

00720 


Xl=FXWAVaP>+FXl5CaPy 

IR 

00730 


X2=FXWAVaPMl>+FX15CUPMl> 

IR 

00740 


FXNET=X1-X2 

IR 

00750 


COOLRaP>=FXNET*FAC/DP 

IR 

00760 

90 

CONTINUE 

IR 

00770 


DO 92 IP=2,NP 

IR 

00780 


IPM1=IP-1 

IR 

00790 


IDD=aPMl*NSB>+l 

IR 

00800 


IDE=IDD-NSB 

IR 

00810 

92 

PRINT 91,PUiaDE> ,PUI aODy ,COOLRaP> 

IR 

00820 

91 

FORMATU0X,2F8.2,4X,F10.2> 

IR 

00830 


STOP 

IR 

00840 


END 

IR 

00850 
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SUBROUTINE INDATA 
C*** READ IN DATA 

DIMENSION GFH20i25,22,ll>,TREi43,15> 

DIMENS ION TAUW ( o , TWIN (1 00 > , TAUU ( 1 00 y , ITX ( 50 > 

DIMENSION TOZON(.43,15,3>,PUU43> 

DIMENSION ISNU140y,FSNai40> 

C0MM0N/H2O/GFH20 

COMMON/PRESS/PUI 

C0Mii0N/C02/TAUW, TAUU, TOZON, TWIN, TRE.tSN 

COiiMON/lNlT/NJPl,NP,NSB,NC,NPT,NW,NT 

NUP1=11 

NP=15 

NPT=43 

NW=22 

NSB=3 

NT=25 

R£AD(5,32>DUMMY 
READ^5,32>DUMMY 
READ0,31> ^PUIUXy ,IX=1,40> 

31 F0RMAT^9F8.3> 

READ ^ 5, 32 > DUMMY 

32 FORMAT (,Al,/,Al,/> 

DO 30 IU=1,NUP1 
DO 30 IW=1,NW 

READO,35y UTXUT) , IT=1 ,NTy 
35 F0RMATU3I6> 

DO 40 IT=1,NT 

GFH20 UT , IW, lU y =ITXUT ) 

40 CONTINUE 
30 CONTINUE 

READi5,32>DUI'iMY 
DO 41 IP=1,NP 

READl,5,42y UTXaX) , IX=1 ,NPT> 

42 F0KMATU4I3y 
DO 43 IX=1,NPT 

43 TREUX,IPy = ITXaX>/10u0. 

41 CONTINUE 

READ C 5, 32 y DUMMY 
READ^5,45yaSNa>,I“1.1140> 

45 F0RMAT140I2) 

DO 50 1=1,1140 
F=ISNUy 

50 FSNa>=“l-E-5*F 
C*** READ TAUW, TAUU 
READ 1,5, 32 y DUMMY 
READi,5,81y CTAUWUWy ,IW=1 ,65y 
READ ^ 5, 32 y DUMMY 
READl5,81y ^AUU UUy , IU=1 , lOOy 
81 FORMATU0F8.5y 
READi5,32yDUMMY 

READl,5,lUyUlT0Z0NaX,IP,LATy ,IP=l,NPy ,IX=l,NPTy ,LAT=l,3y 
10 FORMATU5F5.3y 
READ ^ 5, 32 y DUMMY 
READi5,ly aWINUUy ,IU=l,100y 
1 FORMATUDF4.2y 
RETURN 
END 


IR 00860 
IR 00870 
IR 00880 
IR 00890 
IR 009u0 
IR 00910 
IR 00920 
IR 00930 
IR 00940 
IR 00950 
IR 00960 
IR 00970 
IR 00980 
IR 00990 
IR 01000 
IR OlOlO 
IR 01020 
IR 01030 
IR 01040 
IR 01050 
IR 01060 
IR 01070 
IR 01080 
IR 01090 
IR 01100 
IR OHIO 
IR 01120 
IK 01130 
IR 01140 
IR 01150 
IK 0116u 
IR 01170 
IR 01180 
IR 01190 
IR 01200 
IR 01210 
IR 01220 
IR 01230 
IK 01240 
IR 01250 
IR 01260 
IR 01270 
IR 01280 
IR 01290 
IR 01300 
IK 01310 
IK 01320 
IR 01330 
IR 01340 
IR 0U50 
IR 01360 
IR 01370 
IR 01380 
IR 01390 
IR 01400 
IR 01410 
IR 01420 
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SUBROUTINE GAIxOTO 

IR 

01430 


DIMENSION RRl2,3;,PREU> 

IR 

01440 


DIMENSION NXTaO>,NXBaO;,CSaOj 

IR 

01450 


DIMENSION PUIi43> ,TUIv43; ,PAU43> .RAWU43; ,DPU43> ,TAIl43> 

IR 

01460 


COMMON / CLOUD / CS , NXB , NXT 

IR 

01470 


COMMON/TEMP/TSS 

IR 

01480 


COMMON/PRESS/PUI 

IR 

01490 


COMMON/CLOU/CLTOPllO> ,CLBOTllO> 

IR 

01500 


COMMON/INIT/NUP 1 , NP , NSB , NC , NPT , NW, NT 

IR 

01510 


COMMON/AGIOS/TUI, PAI,TAI,RAWI,UBAR^43j ,WBARl2,43> 

IR 

01520 


C0MM0N/AMHN/TAO5; ,PAi55j ,WA(,55; ,CSSU0> ,NDATA,RLAT,TS,PS 

IR 

01530 


DATA PRE/275.,550./ 

IR 

01540 


DATA RR/I.2164I4,4.62Ill,-0.5I7625E-2,-0.4321E-I,0.1875E-4, 

IR 

01550 


$.11344E-3/ 

IR 

01560 


TSS=TS 

IR 

01570 


DP=lPS-PUH40;>/3. 

IR 

01580 


DO 10 IX=41,43 

IR 

01590 

10 

PUiiix;=PUiax-i>+DP 

IR 

01600 


DO 1010 IX=2,NPT 

IR 

01610 


DPI ax>=tPuiux>-pumx-i>; *1.02 

IR 

01620 


PAiaX>=0.5*lPUIUX^+PJIlIX-l J J 

IR 

01630 

1010 

CONTIi^UE 

IR 

01640 

c 


IR 

01650 

C****K* TEMPERATURE AND HUMIDITY INTERPOLATIONS ************* 

IR 

01660 

c 


IR 

01670 


NCK=*2 

IR 

01680 


DO 1030 IX=2,NPT 

IR 

01690 

1040 

IFv PAIU-> *LT. PAl,NCK>} GO To lU42 

IK 

01700 


IF^NCK .EQ. NDATA> GO TO 1042 

IK 

01710 


NCK=NCK+1 

IR 

01720 


GO TO 1040 

IR 

01730 

1U42 

CONTINUE 

IR 

01740 


XX=ALOG^WAl,NCK-l>; 

IR 

01750 


Xl=ALOG^PAIUX>/PA^NCK-l } ) 

IR 

01760 


X2=ALOGcPUiaX>/PAUJCK-l > / 

IK 

01770 


X3«ALOGl,WAl,NCK.> /WAlNCK-1 ) ) 

IR 

01780 


X4=ALOGlPAlNCK>/PAlNCK-i;> 

IR 

01790 


RAWIUX>=EXP^XX+X1*X3/X4> 

IR 

01800 


IF iRAWiaXy .LT. 0.0> RAWIUX>=U.0 

IR 

01810 


TAiaX>=TAiNCK-U+Xl/X4*iTA^NCK>-TAlNCK-l;; 

IR 

01820 


TUI UX;=TAlNCK-l >+X2/X4*lTAlNCK;-TA(NCK-l ) ) 

IR 

01830 


RAWIUX>=RAWiaX>*DPiaX> 

IR 

01840 

1030 

CONTINUE 

IK 

01850 


TAia> = i.AU2> 

IR 

01860 


TUIU>=TAU2> 

IR 

01870 

C 


IR 

01880 

C**«******* COMPUTE SCALED WATER VAPOR AMOUNTS ************** 

IR 

01890 

C 


IK 

01900 


UBARlU=U. 

IR 

01910 


WBARa,l>=0.0 

IR 

01920 


WBARl2, 1>=0.0 

IR 

01930 


DO 14 IX=2,NPT 

IR 

01940 


DO 12 KB=1,2 

IR 

01950 


SZ 2=RR IKB , 1 )+RR iKB , 2 > *TAI 1 IX>+RR UCB , 3 J *TAI lIX; *TAI ( IX; 

IR 

01960 

12 

WBAR^iCB,IX>=WBARlKB,IX-l;+RAWHIX;*SZ2*^,PAIaX>/PRE^KB;; 

IR 

01970 


XX=iPAHIX;/630.>*RAWiaX;*RAWIiIX;/DPHIXj 

IR 

01980 


XX=XX*EXPUBu0./TAiaX;-6.0811; 

IR 

01990 

14 

UBARUX;=UBARaX-l;+XX 

IR 

02000 


O'" POOR QUAUTY 


c 



IR 

02010 

c********** COMPUTE INDEX LEVELS FOR CLOUDS ****************** 

IR 

02020 

C 



IR 

02030 



DO 20 IC=1,NC 

IR 

02040 


20 

csuc>=cssuc> 

IR 

02050 



DO 3 IC=1,NC 

IR 

02060 



DO 2 IX«=16,42 

IR 

02070 



IF^PUiaX> .LE,CLTOPaC> .AND.PUiaX+ly .GE.CLTOPUC>^GOTO A 

IR 

02080 


5 

IF^PUiaX>.LE.CLBOTUC> .AND.PUIUX+l y .GE.CLBOl UC;>GOTO 6 

IR 

02090 



GOTO 2 

IR 

02100 


4 

UT=CLTOP ( ICy-PUI i IX; 

IR 

02110 



UB=PUIUX+l>-CLTOPaC; 

IR 

02120 



IF^UT.LT.UB>GOTO 7 

IR 

02130 



NXTUC>=IX+1 

IR 

02140 



‘GOTO 5 

IR 

02150 


7 

NXTaC>=IX 

IR 

02160 



GOTO 5 

IR 

02170 


6 

UT=CLBOTaC;-PuiaX> 

IR 

02180 



UB=PUI aX+1 >-CLBO'i ac; 

IR 

02190 



IFiUT.LT.UB>GOTO 8 

IR 

02200 



NXBUG>*=IX+1 

IR 

02210 



GOTO 2 

IR 

02220 


8 

NXBUC;=IX 

IR 

02230 


2 

CONTINUE 

IR 

02240 


3 

CONTINUE 

IR 

02250 



CALL IRH20 

IR 

02260 



CALL RC02031.RLAT; 

IR 

02270 



RETURN 

IR 

02280 



END 

IR 

02290 
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SUBROUTINE IRH20 

DIMENSION IH;43>,DHl43>,SWU; 

DIMENS ION GFH20 1 2 5 , 22 , 1 1> , DELTAl 2 5 ; 

DIMENSION SGU3,15;,SHU5;,SHTU5>,SHSa5> 

DIMENSION RAWlm3; ,TAIl43> ,PAH43> ,TUU43j 
COMMON / AGIOS / TU I , PAI , TAI , RAW I , UBAR I 4 3 > , WEAR ; 2 , 43 J 
C0MM0N/H20/GFH20 
COMMON/TEMP/TS 

COMMON/INIT/NUPl,NP,NSB,NC,NPT,NW,NT 
COMMON / FD I V/ SH , SG , SHT , SHS , FLUXU S 
DATA TEMPI/190. /,SW/-6., -5.4/ 

DATA DW/0.3/, DU/0. 006/, DT/5./ 

DATA DELTA/160., 168., 179., 189., 199., 211., 219., 232., 242., 254., 
$268., 279., 292., 307., 320., 335., 350., 365., 382., 397., 415., 

$432. ,451. ,471. ,492./ 

NU-NUPl-1 

C*** COMPUTE SHS, SHT, SG 
FLUXUS-0. 

DO 43 IP=1,NP 
SHUP7-0. 

SHSaP;=0. 

SHT(,IP>“0. 

DO 44 IX=1,NPT 
SGUX,IP>=0. 

44 CONTINUE 
43 CONTINUE 

FH= iTS-TEMPl > /DT+1 . 501 
■ IG=FH 

IFiIG .GE. NTJ IG=NT-1 

IFlIG .LT. 1^ IG=1 

DS=TS- ITEMP 1+FLOAT^ IG-1 ) *DT) 

NZ=1 

DO 41 KB=1,2 

X=G FH20 lIG , 1 , NZ ;+ VGFH20 iIG+1 , 1 , NZ > -GFH20 UG , 1 , NZ> > /DT*DS , 

FLUXUS-FLUXUS+X 
NZ=NUP1 
41 CONTINUE 

FLUXUS-FLUXUS+DELTAl IGy 
DO 47 IX»1,NPT 
FH=(.TUiaX;-TEMPl>/DT+1.5 
IHUX>=FH 

IF lIHUX; .LT. 1; IHaX;=l 
IF UHUX; .GE. NT> IH(.IX>=NT-1 
F1=IHUX>-1 

DHaX;»TUUIX>-(,TEMPl+Fl*DT> 

47 CONTINUE 

DO 42 IP=1,NP 
IS-aP-l>*NSB+l 
DO 40 IX=1,NPT 
IFUo.EQ.IX>GOTO 40 
IFUX .LT. IS ) IC=IX+1 
IFUX .GT. IS } IC=IX-1 

DY=ABSlOBARaS>-0.5*(UBARaCj+UBARUX>>/ 

KU=DY/DU+1.5 

IF ^KU .GE. NUj KU-NU-1 

FU=DY-FLOAUKU-l >*DU 


IR 

02300 

IK 

02310 

IR 

02320 

IR 

02330 

IR 

02340 

IR 

02350 

IR 

02360 

IR 

02370 

IR 

02380 

IR 

02390 

IR 

02400 

IR 

02410 

IR 

02420 

IR 

02430 

IR 

02440 

IR 

02450 

IR 

02460 

IR 

02470 

IR 

02480 

IR 

02490 

IR 

02 500 

IR 

02510 

IR 

02520 

IR 

02530 

IR 

02540 

IR 

02550 

IR 

02560 

IR 

02570 

IR 

02580 

IR 

02590 

IR 

02600 

IR 

02610 

IR 

02620 

IR 

02630 

IR 

02640 

IR 

02650 

IR 

02660 

IR 

02670 

IR 

02680 

IR 

02690 

IR 

02700 

IR 

02710 

IR 

02720 

Il5. 

02730 

IR 

02740 

IR 

02750 

IR 

02760 

IR 

02770 

IR 

02780 

IR 

02790 

IR 

02800 

IR 

02810 

IR 

02820 

IR 

02830 

IR 

02840 

IR 

02850 
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ORIGINAL Pmt S3 
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C*** KB«1,F0R BAND CENTER REGION, KB*2, FOR BAND WING REGION 
DU 45 KB-1,2 

DX=ABSiWBARUB,lSy-0.5*^WBAR^KB,IC>+WBARlKB,lX>>> 

IF^DX .GT. 30. > GO TO 40 
IF^DX .LT. l.E-6y DX=l.E-6 
DX=ALOG10^DX> 

IFi.DX.LT.SWiKB>>DX=SW^KB> 

KW*^DX-SWtKBy ) /DW+1 . 5001 
IF^^KW .GE. NW> KW=NW-1 
FW=DX-tSWvKB>+FLuATvKW-l y*DW> 

NZ=NUP1 

IFiKB .EQ. 2 ) NZ=KU 
IK=IX 

DO 50 1=1,2 

IF a *EQ- 2 ) IK=IC 

IT=IHaK> 

DHl^DUaKy 

DGDT= ^GFH20UT+1 ,KW,NZ>-GFH20UT,KW,NZ>>/DT 
TDEL=DGDT*DH1 

WDEL=iGFH20UT,KW+l,NZ>-GFH20aT,KW,NZy>/DW*FW 
IF^KB.EQ.2>WDEL=WDEL^-^GFH20aT,KW,NZ+l>-GFH20iIT,KW,NZ> ;/ 
*DJ*FU 

X=GFH20 UT , KW, >+TDEL+WDEL 
IFUK*EQ. l>SHTUP>=SHTaP>+X 

IFaK.EQ.NPT>SHSUP>=SHSaP>+DGDT*^TS-TUI^NPT;> 

IF U-£Q-2> X=*-X 
SGaX,IP>-SGUX,IP>+X 
50 CONTINUE 

45 CONTINUE 
40 CONTINUE 

C 

C******** COMPUTE SHaP> 

C 

IT=lHUi5> 

NZ«1 

DO 46 KB=1,2 

X=GFH20aT,l,NZy + (GFH20aT+l,l,NZ>-GF*i20aT,l,NZ>;/DT*DHaS; 

SHaP>=SHUP>+X 

NZ=NUP1 

46 CONTINUE 

SHaP>=SHaP>+DELTAUT> 

42 CONTINUE 
CALL FDIVER 
RETURN 
END 


IR 

02860 

IR 

02870 

IR 

02880 

IR 

02850 

IR 

02900 

IR 

02910 

IR 

02920 

IR 

02930 

IR 

02 940 

IR 

02950 

Itv 

02960 

IR 

02970 

IR 

02980 

IR 

02990 

IR 

03000 

IR 

03010 

IR 

03020 

IR 

03030 

IR 

03040 

IR 

03050 

IR 

03060 

IR 

03070 

IR 

03080 

IR 

03090 

IR 

03100 

IR 

03110 

IR 

03120 

IR 

03130 

IR 

03140 

IR 

0j150 

IR 

031o0 

IR 

0317U 

IR 

03180 

IR 

03190 

IR 

03200 

IR 

03210 

IR 

03220 

IR 

03230 

IR 

03240 

IR 

03250 

IR 

03260 

IR 

03270 

IR 

03280 

IR 

03290 

IR 

03300 
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OF 


PP;0^ ^ 



SUBROUTINE FDIVER 

IR 

03310 


DIMENSION FLUXUU5^,FLUXDU5>,FUJX^43> 

IR 

03320 


DIMENSION SHTU5> ,SHSU5> ,SHU5> ,SGv43,15) 

IR 

03330 


DIMENSION CSaO),NXBUO;,NXTUO> 

IR 

03340 


COMMON/CLOUD/CS , NXB , NXT 

IR 

03350 


COMMON/INIT/NUP 1 , NP , NSB , NC , NPT , NW, NT 

IR 

03360 


COMMON/ FDIV/ SH , SG , SHT , SHS , FLU XUS 

IR 

03370 


COMMON/ FLUXH/FLUXU , FLUXD 

IR 

03380 


NPM1=NP-1 

IR 

03390 


DO 100 1=1, NP 

IR 

03400 


FLUXUa>“0. 

IR 

03410 


FLUXDa>=0- 

IR 

03420 

100 

CONTINUE 

IR 

03430 

c*** 

COMPUTE DOWNWARD FLUXES 

IR 

03440 


DO 185 IP=2,NP 

IR 

03450 


is=up-i;*nsb+i 

IR 

03460 


ISM1=IS-1 

IR 

03470 


SUM=0- 

IR 

03480 


DO 190 LP=i,ISMl 

IR 

03490 


IX=IS-LP 

IR 

03500 


SUM=SUM+SGUX,IP> 

IR 

03510 


FLUXaX>=SUM 

IR 

03520 

190 

CONTINUE 

IR 

03530 


FLUXa;«FLUXU>-SHTaP>+SHaP> 

IR 

03540 


IFl,NC.NE*0>G0T0 75 

IR 

03550 


FLUXD ap>=FLuxay 

IR 

03560 


GOTO 185 

IR 

03570 

75 

CONTINUE 

IR 

03580 


SUMC=0. 

IR 

03590 


DO 60 IC=1,NC 

IR 

03600 


ncbx=nxbuc> 

IR 

036l0 


IFaS.LE,NXraC>>GOTO 60 

IR 

03620 


IFUS.LE.NCBX>GOTO 62 

IR 

03630 


FLUXD aP>=FLUXDaP>+U’LUXiNCBX;+SHUP>>*CSUC> 

IR 

03640 


GOTO 63 

IR 

U3650 

62 

CONTINUE 

IR 

03660 


FLUXD ap)=FLuxDap>+sHap>*csac> 

IR 

03670 

63 

SUMC=SUMC-K:sac> 

IR 

03680 

60 

CONTINUE 

IR 

03690 


FLUXD UP>=FLUXD UPy+FLUXi U * U • -*SUMC> 

IR 

03700 

185 

CONTINUE 

IR 

03710 

c***< 

COMPUTE UPWARD FLUXES 

IR 

03720 


Du 200 IP=1,NPM1 

IR 

03730 


IS=UP“1>*NSB+1 

IR 

03740 


ISP1=IS+1 

IR 

03750 


SUM=0. 

IR 

0376u 


DO 205 IX=ISP1,NPT 

IR 

03770 


SUM=SUM+SG UX, IP > 

IR 

0j780 


FLUXaX>=SUM 

IR 

03790 

205 

CONTINUE 

IR 

03800 


FLUXiNPT>=FLUX(NPT>+SHSUP>+SHaP> 

IR 

03810 


IFiNC.NE.OlGOTO 76 

IR 

03820 


FLUXUaP>=FLUX^NPT> 

IR 

03830 


GOTO 200 

IR 

03840 

76 

CONTINUE 

IR 

03850 
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SUMC=0. 

IR 

03860 


DO 70 IC*1,NC 

IR 

03870 


NCTX=NXTUC> 

IR 

03880 


IFUS.GE.NXBUC>>GOTO 70 

IR 

03890 


IFaS-GE.NCTX>GOTO 72 

IR 

03900 


FLUXUaP>=FLUXUaP;+iFLUXcNCTX>+SHaP>>*CSaC> 

IR 

03910 


GOTO 73 

IR 

03920 

72 

CONTINUE 

IR 

03930 


FLUXUaP> = FLUXUaP>+SHaP>*CSUC> 

IR 

03940 

73 

SUMC=SUMC+CSaC; 

IR 

03950 

70 

CONTINUE 

IR 

03960 


FLUXUaP>=FLUXUUP>+FLUX^NPT>*U.-SUMC; 

IR 

03970 

200 

CONTINUE 

IR 

03980 


FLUXU^NP>=FLUXUS 

IR 

03990 


RETURN 

IR 

04000 


END 

IR 

04010 
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SUBROUTINE RC0203lRLAT> 

DIMENSION PAI^43^ ,TAIl43> ,RAWI(.43> ,BAI(.44> 

DIMENSION TUU43j,TOZONl43,15,3j 
DIuENS ION TAUU U 00 > , TWIN U 00 > , TAUWl 65 ) 

DIMENSION BWU44; ,BLKWINi25> .BLKC02U5; 

DIMENSION TRE(.43,15/,TAUF(.43,15>,T50JFi43,15> 

DIMENSION FLUXDU5J ,FLUXUa5>,FLUXC«.43> ,FLUXWl43> 

DIMENSION FLUXDCU5> ,FLUXDWU5^ ,FLUXUCU5> .FUIXUWllSy 

DIMENSION CSU0> ,NXBU0> .NXTUOj 

DIMENSION TRU9> ,DTAIU9; ,DTFU9, 19> ,FSNU140> 

COMMON /TEMP/ TS 
COMMON/ FLUXC/ FLU XU , FLU XD 
COMMON/ CLOUD / CS , NXB , NXT 

COMMON / AGIOS/ TU I , PAI , TAI , RAWl , U BAR I 4 3 ; , WuAK ^ 2 > 43 > 
C0Mri0N/C02/TAUW, TAUU , TOZON , TWIN , TRE , FSN 
COMMON/INIT/NUP 1 , NP , NSB , NC , NPT , NW, NT 

DATA TR/266. 6, 262. 82, 257. 13, 251. 35, 245. 54, 239. 95, 234. 54, 
$229.81,227.06,225.19,223.38,221.6,219.8,218.05,216.49, 
$4*216.65/ 

DATA BLKC02/12825. , 14598., 16512., 18568., 20768., 23113., 

1 25603., 28237., 31015., 33937., 37001., 40207., 43552., 

2 47036., 50656., 54412., 58300., 62318., 66465., 70738., 

3 75136., 79655., 84293., 89048., 93918./ 

DATA BLKWIN/ 1944. ,2377. ,2877. ,3450., 4102. ,4838. ,5664. ,6585. , 

* 7606. ,8733. ,9969. , 11320. , 12788. , 14380. , 16097. , 

* 17944. , 19923. ,22038. ,24292. ,26685. ,29221 . ,31902. , 

* 34729. ,37703. ,40825./ 

DATA SW/-4./ ,DW/0. l/,NK/65/,SU/-3.301/,DU/.04/ 

DATA TEMPl/190./,DT/5.0/ 

LOGICAL LINEXP/.TRUE./ 

NPTP1=NPT+1 
LAT=IFIXlRLAT/30.+l . ) 

TAUNPTP1>=TS 
DO 15 IP=1,NPTP1 
DX=TAHIP>-TEMP1 
IT=DX/DT+1 . 5 
IFUT .LT. 2 ) IT=2 
IFUT .GT. NT; IT=NT 
ITM1=IT-1 

X=TEMP l+DT*FLOAT lITMl ; 

DS=TAU1P;-X 

BAI UP 1=BLKC02 I IT;+ I BLKC02 lIT; -BLKC02 UTMl y , *DS/DT 
BWI UP;=BLKW1N UT ;+ {. BLi WIN (.IT; -BLKWIN UTMl ; ; *DS/DT 
15 CONTINUE 

BS=BAI(NPTP1; 

BSS=BWI(,NPTP1; 

C 

C**x****** COMPUTE TRANSMITTANCES IN THE 15 MICRON(TAUF; AND 
C***x***** 9.6 MICRON (TXUF; BANDS. 

C 

DO 25 IP=1,NP 

DO 25 LP=1,NPT 

TXU F(.LP,IP;=TOZON(LP,IP, LAT ; 

25 TAUF(LP,IP;=TREaP,IP; 

IF(. .NOT.LINEXPyGOTO 23 
NSTRA=19 
NSM1=NSTRA-1 
DO 21 LP=1,NSTRA 
21 DTAI(LP;=TAULP;-TR(,LP; 

LL=0 

DO 20 LP=1,NSM1 
LPP1=LP+1 

DO 20 IX=LPP1,NSTRA 
SUh=0. 

IY=IX-LP 


IR 

04020 

IR 

04030 

IR 

04040 

IR 

04050 

IR 

04060 

IR 

04070 

IR 

04080 

IR 

04090 

IR 

04100 

IR 

04110 

IR 

04120 

IR 

04130 

IR 

04140 

IR 

04150 

IR 

04160 

IR 

04170 

IR 

04180 

IR 

04190 

IR 

04200 

IR 

04210 

IR 

04220 

IR 

04230 

IR 

04240 

IR 

04250 

IR 

04260 

IR 

04270 

IR 

04280 

IR 

04290 

IR 

04300 

IR 

04310 

IR 

04320 

IR 

04330 

IR 

04340 

IR 

04350 

IR 

04360 

IR 

04370 

IR 

04380 

IR 

04390 

IR 

04400 

IK 

04410 

IR 

04420 

IR 

04430 

IR 

04440 

IK 

0445u 

IR 

04460 

IR 

04470 

IR 

04480 

IR 

04490 

IR 

04500 

IR 

04510 

IR 

04520 

IR 

04530 

IR 

04540 

IR 

04550 

IK 

04560 

IR 

04570 

IR 

0458o 

IR 

04590 

IR 

04600 

IK 

04610 

IR 

04o2U 

IR 

04630 

IR 

04640 

IR 

04650 

IR 

04660 
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DO 22 L=1,IY 
Ll«LL+1 

SUM=SUM+FSN aL > *DTaI a?+L ) 

22 CONTINUE 
DTFUX,LP>-SUM 
DTFaP,IX>=SUM 

20 CONTINUE 
DO 23 IP=1,7 
IX=aP-l>*NSB+l 
DO 24 LP«1,NSTRA 

24 TAUFaP, IP>=TAUFaP, IP>+DTFaP, IX> 

23 CONTINUE 

DO 28 IP=1,NP 
IS=aP-l>*NSB+l 
DO 32 IX=1,NPT 
IFUS .EQ. IX> GO TO 32 
DX=AB S i WEAR V 2 , IS > -WBAR U • IX > > 
IFKDX.VI. 1 . E-7 >DX=1 . E--7 
IF^DX .LT. 40. > GO TO 35 
X=0. 

YY=0. 

GO TO 40 

35 CONTINUE 
DX=ALOG10^DX> 

IFiDX .LT. SM) GO TO 34 
KW=iDX-SW>/DVH-1.5 
IF(KW .GT. NKy KW=NK 
X=TAUWiKW> 

GO TO 36 
34 X=l. 

36 CONTINUE 

DX=ABS^UBARUS>-UBARaX>> 

IF (DX.LT.l.E-7) GO 10 38 
DX=AL0G10^DX> 

IFiDX .LT. SU> GO TO 38 
KU=lDX-SU>/DU+1.5 
IFiKU .GT. 100> KU=100 
Y=TAUU ^KU } 

YY=TWIN^KU> 

GO TO 40 
38 CONTINUE 
Y=l. 

YY=1. 

40 CONTINUE 

TAUFaX,IP>=»TAUFUX,IP>*X*Y 
TXUF^ IX, IP)=TXUFUX, IP> *YY 
32 CONTINUE 
28 CONTINUE 

Z*** COMPUTE DOWNWARD FLUXES 

fluxjcu>=o. 

FLUXDWa>=0. 

DO 70 IP=2,NP 
FLUXDCaP>*=0. 

FLUXDWUP>=0. 

TAV=1. 

TAYY=1 . 

IS«aP-l>*NSB+l 

ISM1=IS-1 

XC=0. 

XW=0. 


IR 0467U 
IR 046o0 
IR 04690 
IR 04700 
IR 04710 
IR 04720 
IR 04730 
IR 04740 
IR 04750 
IR 04760 
IR 04770 
IR 04780 
IR 04790 
IR 04800 
IR 04810 
IR 04820 
IR 04830 
IR 04840 
IR 04850 
IR 04860 
IR 04870 
IR 04880 
IR 04890 
IR 04900 
IR 04910 
IR 04920 
IR 04930 
IR 04940 
IR 04950 
IR 04960 
IR 04970 
IR 04980 
IR 04990 
IR 05000 
IR 05010 
IR 05020 
IR 05030 
IR 05040 
IR 05050 
IR 05060 
IR Oj)U70 
IR 050 8'J 
IR 05090 
IR 05100 
IR 05110 
IR 05120 
IR 05130 
IR 05140 
IR 05150 
IR 05160 
IR 05170 
IR 0518J 
IR 05190 
IR 05200 
IR 05210 
IR 05220 
IR 05230 
IR 05240 
IR 05250 
IR 05260 
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OF POOR Q'J 




f.3 rrV 

riiwi I ^ 



DO 75 LP=1,ISM1 

IR 

05270 


IX«IS-LP 

IR 

05280 


XC»XC+BAiaX+l ;*vTAY-TAUFUX, IF) j 

IR 

05290 


XW=XW+BWI K IX+1 ) * ( TAYY-TXU F UX , IP > > 

IR 

05300 


FLUXCUX>«XC 

IR 

05310 


FLUXWaX;=XW 

IR 

0532o 


TAY«TAUFUX,Ii:'> 

IR 

05J30 


TA Y=TXUFUX,IP> 

IR 

05340 

75 

CONTINUE 

IR 

05350 


IF^NC.NE.O>GOTO 69 

IR 

05360 


FLU XDC a P > “FLc XC a > 

IR 

05370 


FLUXDWUP>=FUJXW^ 1 ) 

IR 

05380 


GOTO 70 

IR 

05390 

69 

CONTINUE 

IR 

05400 


SUMC=0. 

IR 

05410 


DO 60 LC=1,NC 

IR 

05420 


NCBX=NXBUC> 

IR 

05430 


IFaS.LE.NXraC>)GOrvJ 60 

IR 

0544U 


IFUS.LE.NCBX>GOTO 62 

IR 

05450 


FLUXDCaP>=FLUXDCUF>+^FLUXC^NCBX>+BAI^NCBXy*TAUF^NCBX, 

IPyy*CSaCyIR 

05460 


FLUXDWUP> = FLUXDWUP>+^FLUXWi,NCBXy+BWl^NCBX>*TXUF^NCBX, 

IPyy*CSUCyIR 

05470 


GOTO 63 

IR 

05480 

62 

CONTINUE 

IR 

05490 


FLUXDCUP>=rUJXDCUP>+BAiaS^*CSaC> 

IR 

05500 


FLUXDWaP>=FLUXDWaP>+BWiaSy*CSUCy 

IR 

05510 

63 

SUMC=SUMC+CSUC> 

IR 

05j20 

60 

CONTINUE 

IR 

05530 


FLUXDClIPy=FLUXDClIPy+FLUXCUy*ll*-SUMCy 

IR 

05540 


FLUXUWaPy=FLUXDWaP>+FLUXWa >*i 1 •-SUMC; 

IR 

05550 

70 

CONTINUE 

IR 

05560 

c*** 

COMPUTE UPWARD FLUXES 

IR 

0557U 


NPMi=NP-l 

IR 

05580 


DO 80 IP=1,NPM1 

IR 

05590 


FLUXUCUP>=0. 

IR 

05600 


FLUXUWUP>=0. 

IR 

05610 


TAY=1 . 

IR 

05620 


TAYY=1. 

IR 

05630 


IS=aP'ly*NSB+l 

IR 

05640 


ISP1=IS+1 

IR 

05650 


XC=0. 

IR 

05660 


xw=o. 

IR 

05670 


DO c5 IX=ISP1,NPT 

IR 

05680 


XC=XC+BAI UX> * TAY-TAUF UX , IP y > 

IK 

05690 


XW=XW+BW I a X y * TAYY-TXU F UX , IP y y 

IR 

05700 


FLUXCUXy=XC 

IR 

05710 


FLUXWiIX>=XW 

IR 

05720 


TAY=TAUFUX,IPy 

IR 

05730 


TAYY=TXUFUX,lPy 

IR 

05740 

85 

CONTINUE 

IR 

05750 


FLUXCiNPTy=FLUXC^NPTy+BS*TAY 

IR 

05760 


FLUXW^NPTy=FLUXW^NPT>+BSS*TAYY 

IR 

05770 


IF^NC.NE.OyGOTO 68 

IR 

05780 


FLUXU C ( IPy =FLUXC 1 NPT y 

IR 

05790 


FLUXUWUP>=FLUXW^NPTy 

IR 

05800 


GOTO 30 

IR 

05810 

66 

CONTINUE 

IK 

05620 


SUMC=0, 

IR 

05830 


DO 64 IC=1,NC 

IR 

05840 


NCTX=Nxrac> * 

IR 

05850 


IFUS.GE.NXBUCyyGOTO 64 

IR 

05860 


IFUS .Gh . NCTXyGOTO 66 

IR 

05870 


A-12 


OF POOR 


QUAirry 


FLUXUC{lPj°FUJXJC{lP)+{FWXC(SCTX)+BAl{UCTX+l)*TAUF{NCTX,lPj) 

i*csac> 

FLUXUWUP>“FLUXUWaP>+vFLUXW(NCTX>+BWIvNCTX+l>*TXUF^NCTX,IP> 

1>*CSUC) 

GOTO 67 

66 CONTINUE 

FLUXUCaP>=FLUXUCaP>+BAiaS+l>*CSUC> 
FLUXUWaP>“FLUXUWaP>+BWIUS+l >*CSUC> 

67 SUMC=SUMC+CSUC> 

64 CONTINUE 

FLUXUCUP>=FLU:ajCaP>+FLUXC^NPT>*a*^SUMC> 

FLUXUWUP>“FLUXUWaP>+FLUXW^NPT>*U*-SUMC> 

80 CONTINUE 

FLUXUC^NP>«BS 

FLUXUW^NP>«BSS 

C 

C********* total fluxes in the 9.6 AND 15 MICRON BAt^S ************* 
C 

Do 97 1P=1,NP 

FLUXU aP>“FLUXUCUP>+FLUXUW^ IP> 

FLUXD aP>=FLUXDCc IP>+FLUXDWi IP ) 

97 CONTINUE 
RETURN 
END 


IR 05880 
IR 05890 
IR 05900 
IR 05910 
IR 05920 
IR 05930 
IR 0594U 
IR 05950 
IR 05960 
IR 05970 
IR 05980 
IR 05990 
IR 06uu0 
IR 06010 
IR 06020 
IR 06030 
IR 06040 
IR 06050 
IR 06060 
IR 06070 
IR 06080 
IR 06090 
IR 06100 
IR UoliO 
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APPENDIX B 


ARRAY AND VARIABLE LIST 


I . Main Program 


Code 

Meaning 

Unit 

PA 

Input pressure level 

mb 

TA 

Input temperature 

K 

WA 

Input specific humidity 

gm/gm 

CLTOP 

Cloud top pressure 

mb 

CLBOT 

Cloud base pressure 

mb 

CSS 

Fractional cloud amount 


NDATA 

Number of input data in PA, TA, WA 


NC 

Number of cloud types 


RLAT 

Latitude 

degree 

PS 

Surface pressure 

mb 

TS 

Surface temperature 

K 

PUI 

Pressure sublevel 

mb 

FLUXHD 

Downward flux in the H 2 O bands 

ergs/(cm2 sec) 

FLUXHU 

Upward flux in the H 2 O bands 

ergs/(cm2 sec) 

FLUXCD 

Downward flux in the 9.6 and ISpm bands 

ergs/(cm2 sec) 

FLUXCU 

Upward flux in the 9.6 and 15/am bands 

ergs/(cm2 sec) 

FXWAV 

Net flux in the H 2 O bands 

ergs/(cm^ sec) 

FX15C 

Net flux in the 9.6 and 15pm bands 

ergs/(cm^ sec) 

FAC 

Conversion factor for cooling rate 


COOLR 

Cooling rate 

°C/day 
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2. Subroutine INDATA 


Code 

Meaning 

Unit 

PUI 

See main program 


GFH20 

G-function 

ergs/(cm2 sec) 

TRE 

CO 2 transmission function 


FSN 

Coefficients for correcting CO 2 transmittance 

K-i 

TAUW 

Transmittance in the ISjum band due to H 2 O line absorption 


TAUU 

Transmittance in the 1 S^rm band due to H 2 O e-type absorption 


T020N 

Ozone transmission function 


TWIN 

Transmittance in the 9.6/im band due to e-type absorption 


NP 

Number of flux levels 

1 

NPT 

Total number of sublevels 


NSB 

Number of sublayers between neighboring flux levels 


NT 

Number of temperature values in GFH20 


NW 

Number of log w values in GFH20 


NUPl 

Number of u values in GHF20 
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3. Subroutine GAMOTO 


Code 

Meaning 

Unit 

PA 

See main program 


TA 

See main program 


WA 

See main program 


PUI 

See main program 


PAI 

Pressure at the middle of sublayers 

mb 

DPI 

Air mass in a sublayer 

gm/cm^ 

TUI 

Temperature at sublevels 

K 

TAI 

1 

Temperature at the middle of sublayers 

K 

RAWI 

Specific humidity in a sublayer (01800) 

gm/gm 

RAWI 

Water vapor content in a sublayer (01840) 

gm/cm^ 

TS, TSS 

Surface temperature 

K 

WEAR 

Scaled water vapor amount for line absorption 

gm/cm^ 

UBAR 

Scaled water vapor amount for e-type absorption 

gm/cm^ 

PRE 

Reference pressure 

mb 

RR 

Coefficients for computing R 


CS, CSS 

Fractional cloud cover 


NXT 

Index for cloud top pressure 


NXB 

Index for cloud base pressure 


CLTOP 

See main program 


CLBOT 

See main program 

1 

NC 

See main program 
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4. Subroutine IRH20 


Code 

Meaning 

Unit 

TEMPI 

Initial temperature value for array GFH20 

K 

SW 

Initial value of log w for array GFH20 


DT 

Increment in temperature values in GFH20 

K 

DW 

Increment in log w in GFH20 


DU 

Increment in u in GFH20 

gm/cm^ 

IH 

Temperature index for GFH20 


KW 

Index of log w for GFH20 


KU 

Index of u for GFH20 


DH 

Difference between the sublevel temperature and the temperature 
at the nearest index value. 

K 

FW 

Difference between the value of log w and that at the nearest 
index value. 


FU 

1 

Difference between the value of u and that at the nearest index 
value. 

gm/cm^ 

SH 

Planck flux at sublevels 

j 

ergs/(cm2 sec) 

FLUXUS 

Upward flux at the surface 

ergs/(cm2 sec) 

SHT 

G(p, pt, Tj) 

ergs/(cm2 sec) 

SHS 

G(p, Ps, Tj) - G(p, pj, T(ps)) 

ergs/(cm^ sec) 

SG 

OG/aT) AT 

eigs/(cm2 sec) 

GFH20 

See subroutine INDATA 


DELTA 

Correction term for Planck fluxes 

ergs/(cm^ sec) 

DGDT 

dG/8T 

ergs/(cm2 secK) 

TDEL 

OG/9T) AT 

ergs/(cm^ sec) 

WDEL 

(9G/9 log w) A log w, (9G/9u) Au 

ergs/(cm2 sec) 

NUPl 

See subroutine INDATA 


NP 

See subroutine INDATA 



B-4 





4. Subroutine IRH20 (Continued) 


Code 

Meaning 

Unit 

NPT 

See subroutine INDATA 


NSB 

See subroutine INDATA 


NW 

See subroutine INDATA 


NT 

See subroutine INDATA 


TS 

See main program 


TUI 

See subroutine GAMOTO 


WBAR 

See subroutine GAMOTO 


UBAR 

See subroutine GAMOTO 
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5. Subroutine FDIVER 


Code 

Meaning 

Unit 

FLUX 

2 SG 

ergs/(cm2 sec) 

FLUXU 

Upward flux in the water vapor bands 

ergs/(cm2 sec) 

FLUXD 

Downward flux in the water vapor bands 

ergs/(cm2 sec) 

SUMC 

Total cloud cover 


NP 

See subroutine INDATA 


NPT 

See subroutine INDATA 


NSB 

See subroutine INDATA 


NC 

See main program 


cs 

See subroutine GAMOTO 


NXT 

See subroutine GAMOTO 


NXB 

See subroutine GAMOTO 


SG 

See subroutine GAMOTO 


SH 

See subroutine GAMOTO 


SHS 

See subroutine GAMOTO 


SHT 

See subroutine GAMOTO 


FLUXUS 

See subroutine GAMOTO 
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6. Subroutine RC0203 


Code 

Meaning 

Unit 

TEMPI 

Initial temperature value for arrays BLKC02 and BLKWIN 

K 

SW 

Initial value of log w for TAUW 


SU 

Initial value of log u for TAUU and TWIN 


DT 

Increment in temperature values in BLKC02 and BLKWIN 

K 

DW 

Increment in log w in TAUW 


DU 

Increment in log u in TAUU and TWIN 


NT 

Number of temperature values in BLKC02 and BLKWIN 


NK 

Number of log w values in TAUW 


LINEXP 

Flag for correcting CO 2 transmittance due to temperature 
variations 


RLAT 

Latitude 

degree 

LAT 

Index for latitude 


TAl 

See subroutine GAMOTO 


TR 

U.S. Standard Temperature profile 

K 

DTAI 

Temperature deviation from TR 

K 

WEAR 

See subroutine GAMOTO 


UBAR 

See subroutine GAMOTO 


BLKC02 

Precomputed Planck flux in the 1 Spm band 

ergs/(cm2 sec) 

BLKWIN 

Precomputed Planck flux in the 9.6pm band 

ergs/(cm^ sec) 

BAI 

Planck flux in the middle of sublayers m the 15pm band 

ergs/(cm2 sec) 

BWI 

Planck flux in the middle of sublayers in the 9.6pm band 

i 

ergs/(cm2 sec) 

BS 

Upward flux at the surface in the 15pm band 

ergs/(cm2 sec) 

BSS 

Upward flux at the surface in the 9.6pm band 

ergs/(cm2 sec) 

ESN 

See subroutine INDATA 


TOZON 

See subroutine INDATA 
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6. Subroutine RC0203 (Continued) 


Code 

Meaning 

Unit 

TRE 

See subroutine INDATA 


TAUW 

See subroutine INDATA 


TAUU 

See subroutine INDATA 


TWIN 

See subroutine INDATA 


TAUF 

Transmittance in the ISjum band 


TXUF 

Transmittance in the 9.6fim band 


FLUXC 

Z B A T in the 1 5/im band 

ergs/(cm2 sec) 

FLUXW 

2 B A T in the 9.6/^m band 

ergs/(cm2 sec) 

FLUXDC 

Downward flux in the 1 5/im band 

ergs/(cm2 sec) 

FLUXUC 

Upward flux in the 1 5pm band 

ergs/(cm2 sec) 

FLUXDW 

Downward flux in the 9.6pm band 

ergs/(cm^ sec) 

FLUXUW 

Upward flux in the 9.6pm band 

ergs/(cm2 sec) 

FLUXU 

Total upward flux in the 9.6 and 15pm bands 

ergs/(cm2 sec) 

FLUXD 

Total downward flux in the 9.6 and 15pm bands 

ergs/(cm2 sec) 

NC 

See main program 


NXT 

See subroutine GAMOTO 


NXB 

See subroutine GAMOTO 


CS 

See subroutine GAMOTO 


SUMC 

Total cloud cover 


NP 

See subroutine INDATA 


NPT 

See subroutine INDATA 


NSB 

See subroutine INDATA 
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FAGu u; 

OF POOR QUALITY 


APPENDIX C 
SAMPLE OUTPUT 

1. Clear Mid-Latitude Winter Atmosphere (NC = 0) 


NUATA NC KLAT TS PS 

0 15.0 270.0 1013.0 


PRESSURE LEVELS (.MB 

j OF SAMPLE 

PROFILE 


0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.5 

10.0 

12.5 

15. U 

17.5 

2u.0 

25.0 

30.0 

35.0 

40.0 

50.0 

60.0 

70.0 

85.0 

100.0 

125.0 

150.0 

175.0 

200.0 

250.0 

300.0 

350.0 

400.0 

500.0 

600.0 

700.0 

850.0 

920.0 

1000.0 

IFMPERATURESlK/AT 

THE ABOVE PRESSURE LEVELS 


240.64 

252.60 

258.73 

264.03 

263.21 

255.07 

243.47 

236.23 

230.82 

226.98 

224. 7u 

222.30 

219.89 

217.36 

216.80 

216.33 

215.93 

215.42 

215.20 

215.20 

215.20 

215.20 

215.20 

215.31 

215.87 

216.43 

217.05 

217.69 

218.23 

218.69 

219.27 

222.39 

229.03 

234.55 

242.40 

251.24 

258.73 

264.69 

268.32 

270.57 

HUMIDITY IGM/GM/ Ai 

THE ABOVE 

PRESSURE LEVELS 


0.285E-05 

0.439E-05 

0.548L-05 

0.663L.-05 

0.820L-0^ 

O.lllE-04 

0. 118E-04 

0.127E-04 

0.135E-04 

0. 143E-04 

0.135E-04 

0. 169E-04 

0.185E-04 

0.202E-04 

0. 193C.-04 

0.18bE-04 

0.180E-04 

0. 173E-04 

0.138E-04 

0- 105E-04 

0.853E-05 

0.700E-05 

0.512E-05 

0.4O8E-05 

0.415E-05 

0.400L-05 

O.400E-05 

0.453E-05 

0.745E-05 

0.195E-04 

0.195E-04 

U.244E-04 

0.500E-04 

0. lOUE-03 

0.272E-03 

0.574L-03 

0.102E-02 

0.166E-02 

0.210E-02 

0.242E-02 


CLOUD TOP CLOUD BOTTOM CLOUD FRAuTION 
210. 300. 0. 18 


PRESSURE FLU XL:. IN H20 BANDS FLuXES IN 9- 6 & 15 MIv^RON BANDS 


iMby ^ERGS/CM**2/SEC> ^ ERGS/CM**2/S£Cj 



OuWN 

UP 

DOWN 

UP 

1.39 

0.0 

185444.37 

0.0 

48202.96 

3.18 

267.87 

185j9u.3i 

2081.45 

47522.17 

7.28 

287.19 

185349.50 

2958.25 

46955.96 

16.67 

632.25 

185324.37 

4343.21 

47013.10 

38.16 

1855.69 

185351.94 

t)837.71 

48156.27 

8/. 36 

2985.56 

185424.69 

10492.16 

500/2.78 

200.00 

6850.19 

185616.37 

14897.36 

5197:,. 02 

300.00 

14780.50 

187467.50 

17737.02 

54861.20 

4UU.00 

28260.62 

190267.00 

23370.52 

59365.65 

500.00 

47726.87 

194613.50 

29131.07 

63406.77 

600.00 

69482.19 

20U69.J6 

35580.37 

67895.44 

700.00 

92101.81 

208026.37 

42254.19 

71718.94 

8u0.00 

110372.37 

213865.00 

47685.18 

74349.87 

900.00 

12806U.00 

219744.00 

52466.74 

77-83.87 

1013.00 

147111.94 

223182.00 

57425.71 

78223.00 


PRESSURE 

COOLING RATE 



lMBAR> 

^DEGREE CLLCIUS/DAY> 



FROM 

TO 



1.39 

3 

.18 

14. 

.55 

3.18 

7 

.28 

3. 

.09 

7.28 

16 

.67 

1. 

.71 

16.67 

38 

. 16 

0. 

.92 

38.16 

87 

.36 

0. 

.48 

87.36 

200 

.00 

0. 

. 46 

200.00 

300 

.00 

0. 

.51 

300. 00 

400 

.00 

1. 

.00 

400.00 

500 

.00 

1. 

.42 

SoO.OO 

600 

.00 

1. 

.44 

600.00 

700 

.00 

1. 

.58 

700.00 

ooO 

.00 

1. 

.29 

SoO.OO 

900 

.00 

1. 

.16 

900.00 

1013 

.00 

1. 

.46 
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ORIGR>JAL 

OF POOR QUALirf 


2. Cloudy Mid-Latitude Winter Atmosphere (NC = 5) 



NDATA NC RLAT 

TS PS 




40 5 15.0 

270.0 lul3 

.0 



PRESSURE LEVELalMBy 

OF SAMPLE 

PROFILE 


0-2 

0.3 

0.5 

0.7 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.U 

8.5 

10.0 

12.5 

15.0 

17.5 

20.0 

25.0 

30.0 

35.0 

40.0 

50.0 

60.0 

70.0 

85.0 

luO.u 

125.0 

150.0 

175.0 

200.0 

250.0 

300.0 

350.0 

400.0 

500.0 

60U.0 

700-0 

850.0 

920.0 

1000.0 


TEMPERATURES iK>AT THE ABOVE PRESSURE LEVELS 


240.64 

252.60 

258.73 

264.03 

263.21 

255.07 

243.47 

236.23 

230.82 

226.98 

224.70 

222.30 

219.89 

217.36 

216.80 

216.33 

215.93 

215.42 

215.20 

215.20 

215.20 

215.20 

215.20 

215.31 

215.87 

216.43 

217.05 

217.69 

218.^3 

218.69 

219.27 

222.39 

229.03 

234.55 

242.40 

251.24 

258.73 

26^.69 

266.32 

270.57 

HUMIDIlYlGi./GMy 

AT THE ABOVE 

PRESSURE LEVELS 


0.285t-05 

0.439E-05 

0.548E-05 

0.663E-05 

0.820E-05 

0. 11 IE-04 

0.118E-O4 

0.127E-04 

0.135E-U4 

0. 143E-04 

0.155E-04 

0.169E-04 

0.185E-04 

0.202E-04 

0.193L-04 

0.186E-04 

0.180E-04 

0.173E-04 

0.138E-04 

0.105E-04 

0.C53E-05 

0.700E-05 

0.512E-05 

0.468E-05 

0.415E-05 

0.400E-U5 

O.mOOL-05 

0.45JE-05 

0.745E-05 

0.195E-04 

0.195E-04 

0.244E-04 

0.500E-04 

O.lOOE-03 

0.272E-03 

0.574E-03 

0.102E-02 

0.166E-O2 

0.210E-02 

0.24 2 £-02 


CLOUU TOP CLOUD BOTTOM CLOUD FRACTION 


210. 

300. 

0. 18 

5u5. 

600. 

0.06 

601. 

901. 

0. 15 

6uU. 

900. 

0.06 

204. 

902. 

0.04 


PuESbURE 

FLUXES 

IN H20 BANDS 

FLUXES IN 9. 6 

6 15 MICRON BANDS 

IMB> 

URGS/CM**2/SEC> 

URGS/CM* 

*2/SEC> 


DOWN 

UP 

DOWN 

UP 

1.39 

0.0 

160269.50 

u.o 

43336.52 

3.18 

267.87 

160215.44 

2081.45 

42645.90 

7.28 

287.19 

16017^.56 

2958.25 

42049.80 

16.67 

832.25 

160148.75 

4343.21 

41992.21 

38.16 

1855.69 

160173.62 

6o37. 71 

42797.15 

87.36 

2985.56 

160240.19 

10492.16 

44245.09 

2uO.00 

6550. 19 

160410.25 

148‘^7.36 

45780.01 

3 Jo. 00 

34820.18 

177958.44 

20919.61 

52523.17 

400.00 

46781-64 

181566.06 

26296.43 

57142.21 

500.00 

63775.48 

186278.69 

31799.63 

61207.65 

600.00 

89371.12 

196293.62 

39297.41 

66551.94 

70U.0O 

Il5u26.87 

204333.69 

46916. '.8 

70682.81 

800.00 

130624.25 

210515.56 

51979.39 

73420.44 

9o0.00 

159386.06 

219744.00 

59768.91 

77283.87 

013.00 

173253. j7 

223182.00 

6^878.33 

78223.00 


PRESSURE COOLINC RATE 

KKBARj ^DEGREE CLLCIUS/DAY> 

FROM TO 


1. 

.39 

3. 

.18 

14. 

.59 

3. 

.18 

7. 

.28 

3. 

.16 

7. 

.28 

16. 

.67 

1. 

.81 

16. 

.67 

j8. 

.16 

1. 

.06 

38. 

.16 

87. 

.36 

0. 

.56 

87, 

.36 

200- 

.00 

0. 

.49 

20u, 

.00 

30o. 

.00 

0, 

.82 

300. 

.00 

400. 

.00 

0. 

.77 

400. 

.00 

500. 

.00 

1. 

.16 

500, 

.00 

6oO, 

.00 

1. 

,5o 

600. 

.00 

700, 

.00 

1. 

.78 

700, 

OU 

800, 

.00 

0, 

.99 

800, 

.00 

900. 

.^0 

I, 

.98 

90u. 

Ou 

1013. 

00 

I, 

.02 
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ORIGIPiAi PASS Ki 

OF POOR QUALITY 

3. Clear Tropical Atmosphere (NC = 0) 

NDATA NC RLAI TS PS 

4U 0 15.0 297.8 1013.1 

PRESSURE LEVELS OF SAMPLE PROFILE 


0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

;.u 

b.5 

10.0 

12.5 

15.0 

17.5 

20.0 

25.0 

30.0 

35.0 

40.0 

50.0 

60.0 

70.0 

85.0 

100.0 

125.0 

150.0 

175.0 

200.0 

25 n.u 

300.0 

350.0 

400.0 

5UO.O 

600.0 

700.0 

850.0 

920.0 

1000.0 


TEMPERATURES UMT 

THE ABOVE 

PRESSURE LEVELS 


229.36 

246.72 

255.63 

263.31 

269.91 

265.47 

257.36 

251.76 

247.68 

244.41 

241.76 

239.03 

236.29 

233.26 

230.23 

227.77 

225.66 

222.96 

220.12 

217.96 

216.11 

212.61 

<^07.6c 

203.59 

199.42 

195.30 

197.54 

205.47 

211.85 

218.32 

22 b . 22 

235.13 

243.39 

249.66 

259.20 

269.24 

278.52 

286.72 

292.90 

297.17 


HUMIDITYiGM/GM> AT 

THE ABOVE 

PRESSURE LEVELS 


0.200t-U5 

0.3l4h-05 

0.395E-U5 

0.482E-05 

0.570E-05 

0.870E-05 

0.990E-05 

0.1U7E-04 

0.117E-04 

0. 124E-04 

0. 134E-04 

0.149E-04 

0.166E-04 

0. 187E-04 

0.191E-04 

0.184E-U4 

u. 178I.-04 

0.171E-04 

0.147E-04 

0.109E-U4 

0.871E-05 

0.692E-05 

0.498E-05 

0.435E-05 

0.368E-05 

0.336E-05 

0.327E-05 

0.349E-05 

0.437E-U5 

0.755E-05 

0.255E-04 

0.919E-U4 

0.239E-03 

0.458E-03 

0.936E-03 

0. 196E-02 

0.319E-02 

0.794E-02 

0. 116E-01 

0.141E-01 


CLOUD TOP CLOUD BOTTOM CLOUD FRACTION 
200. 300. 0.45 


PRESSURE FLUXES IN H20 BANDS FLUXES IN 9- 6 & 15 MICRON BANDS 


^MB^ URGS/CM**J/SEC> URGS/Cm**2/SEC; 



DOWN 

UP 

DOWN 

UP 

1.39 

0.0 

236907.69 

0.0 

62255. 7j 

3.18 

30d.94 

236876.50 

2649.86 

61631.22 

7.28 

307.37 

2368^7.25 

4193.17 

60845.29 

16. 67 

933.00 

236673.75 

6140.75 

60214.68 

38.16 

2034.50 

236570.31 

8110.90 

607j0.78 

87.36 

2685.94 

236504.87 

8774.48 

62061.51 

200.00 

5388.50 
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